Static and dynamic image potential for tunneling into a Luttinger liquid 
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We study electron tunneling from a tip or a lead into an interacting quantum wire described by 
Luttinger liquid theory. Within a WKB-type approach, the Coulomb interaction between the wire 
and the tunneling electrons, as well as the finite traversal time are taken into account. Although the 
static image potential is only logarithmically suppressed against the bare Coulomb interaction, the 
dynamic image potential is not strong enough to alter power-law exponents entering the tunneling 
density of states. 



PACS numbers; 7f.iO.Pm, 72.80.Rj, 73.40. Gk 

One-dimensional (fD) quantum wires (QWs) are at 
the focal point of current activities in condensed- 
matter physics. pFabrication advances iit, semiconductor 
heterostructuresEl and carbon nanotubesci allow the sys- 
tematic study of phenomena arising only in one dimen- 
sion. In particular, the physics of ID nanowires is in- 
timately connected with the concept of a Luttinger liq- 
uid (LL), exhibiting spin-charge separation, suppression 
of the tunneling density of states (TDOS), and inteij. 
action dependent power-laws in transport properties E'u 
For nanotubes, the theoretically predictedEl LL behav- 
ior has been-, convincingly established in several recent 
experiments.l3 Of particular importance was the measure- 
ment of the TDOS power-law exponent a, which provides 
information about the dimensionless interaction strength 
parameter g of the LL (where < 5 < 1). The TDOS 
can be obtained, e.g. from the tunneling current mea- 
sured via weakly attached metallic leads. Another pos- 
sibility is to use a scanning tunneling microscopy (STM) 
tip. Apart from the resolution of atomic and electronic 
properties of individual nanptubes successfully achieved 
in recent STM experiments,Lru a detailed analysis of non- 
linear current-voltage curves in the spectroscopy (STS) 
mode would allow to extract the TDOS exponent. 

So far in all theoretical studies of tunneling into a LL, 
it has tacitly been assumed that the electrons tunnel- 
ing into the correlated fermion liquid do laojb modify the 
effective TDOS. Conventional treatmentaij employ the 
tunneling Hamiltaman and thereby assume the traversal 
time of tunnelingHil to be infinitely short. In addition, 
the interaction between the tunneling electron and the 
LL electrons is neglected. Given these two assumptions, 
the TDOS is indeed an intrinsic property of the LL with 
a power-law dependence on energy, p{e) (x e", for e ^ 0, 
so that a depends only on g. (For clarity, we focus on 
the case of bulk tunneling at zero temperature through- 
out this paper.) For an STS experiment, such a calcula- 
tion predicts that the measured differential conductance 
is proportional to the TD0S,t3 and hence the exponent 
a can be extracted from experimental data. A similar 
reasoning has been employed to understand the data in 
Ref. g where electron tunneling between metallic leads 
and a nanotube was important. 



Here we address the question of whether this measured 
exponent reaUy characterizes the unperturbed LL, or 
whether it is affected by the dynamics and charge of the 
tunneling particle. To that end, we take into account cor- 
relations both within the QW and between the QW and 
the tunneling electron, as well as a finite traversal time. 
Our main findings are as follows: Although the static 
image potential experienced by the tunneling electron is 
very strong, dynamical effects turn out to be of crucial 
importance. Within the framework of a semiclassiGal the- 
ory related to Nazarov's tunnel junction theory,Ej the 
low-energy power-law exponent of the TDOS is governed 
solely by intrinsic LL properties. Thereby we provide 
the a posteriori justification for the (previously assumed) 
connection linking the (observable) value of a and the in- 
teraction parameter g. Since spin and charge are decou- 
pled in a LL, we study only the spinless single-channel 
case, with the same conclusions applying to spin-i elec- 
trons or nanotubes. 

For the case of tunneling from an STM tip into the 
QW, the relevant geometry is depicted in Figure |l|. We 
consider a clean and very long (L —f 00) QW, and mainly 
focus on an interaction potential of the form 
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v/x2 + {D- vY + a2 



(1) 



where a is the lattice constant of the QW. For y — D, 
this describes the intra-wire interaction responsible for 
the LL state, while for < y < -D, it gives the interaction 
between the wire electrons and the tunneling electron at 
X = (see below). The one-sided Fourier transform of 
Eq. d) is 



Viq,y)^2Ko{\qWiD-yr+a' 



(2) 



with the modified Bessel function Kq. When using 
Eq. (^ , one neglects or crudely approximates several im- 
portant effects, e.g. screening due to the tip, or the orbital 
structure of both tip and quantum wire (see, e.g. Ref. |^). 
However, our qualitative findings concerning the connec- 
tion between a and g are expected to hold quite generally. 
In addition to the model potential (m) , these findings can 
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be made rigorous for a general class of separable poten- 
tials. 

Tunneling out of the tip (or lead) proceeds from the 
ground state with energy E'q into the LL ground state 
with energy Eq, and is usually described in terms of a 
simple tunneling Hamiltonian, 



Ht = r«'t(0)co + h.c. 



(3) 



where cq annihilates an electron in a state at the center of 
the tip and ^'^^(0) creates an electron at a; = in the LL. 
If the tunneling matrix element T is small, Ht can be 
treated as a perturbation to H' + H, where H' describes 
the electrons in the tip and H denotes the Hamiltonian of 
the LL. More specifically, the perturbative treatment is 
appropriate when the tunneling resistance is large com- 
pared to h/e^. Then the tunneling rate F can be calcu- 
lated with the aid of Fermi's golden rule. Labeling states 
in the tip by \ v) and states in the LL by \n), the tunneling 
rate is (we put h = 1) 



27rT2^|(z.,n|vI/t(0)co|0,0>p<5(i?^ 



-Eq- E'^- En) , 



where |0, 0) denotes the ground state in the absence of 
tunneling. Using the identity 

SiE'a+Eo-El-E„) ^ j de 5{Eo+t-En)5{E'o-e-El) , 



the rate may be written as 

sV 



T = 2tiT^ I de p{e)p'{eV - e) 
Jo 



(4) 



Here the TDOS for adding an electron with energy fi + e 
to the LL is 

p{e) = ^ |(n|vl.t(0)|0)|25(i5„ +p + e- E,,) , (5) 

n 

where the electrochemical potential fi is the minimal en- 
ergy required to add an electron. Furthermore, 

p'(e') = ^ |(Hco|0)p5(i?^ -p' + e'- El) (6) 



is the DOS for removing an electron from the tip/lead 
with energy p' — e' . Now p{e) and p'{e) are nonvanishing 
only for e > 0, and p' — p = eV determines the ap- 
plied voltage. With these definitions we readily obtain 
Eq. (||). In a metallic lead, the DOS p'{e) is essentially 
constant, while for an STM tip, we expect pronounced 
peaks reflecting the discrete level structure. In the latter 
case, Eq. (||) reproduces the rate from an STM^tip into 
a metal first obtained by Tersoff and Hamann.lla 
To determine p(e) explicitly, we write 



p{e) = — dt y(0|*(0)|n)(n|*t(o)|0)e*(^'' 



Re 



dt Go(t)e'(^+')* , 



(7) 



with the single-electron Greens function (t > 0) for an 
electron at position a; = in the ground state, G'o(t) = 
(^'(0,t)^''l'(0,Qi^. To evaluate Go(i), standard bosoniza- 
tion methodsEfQ can be applied. The kinetic part of the 
Hamiltonian is 



+ (d^if)^] + pN 



Here is the Fermi velocity, N the particle number 
operator, and 19(2;) is conjugate to the phase field f{x) 
describing the plasmon excitations in the wire. The 
intra-wire interaction part is Hy = i Jdxdx' n{x)V{x~ 
x' , D)n{x') with the electron charge density n(x) = 
—eTr~^^'^dx'P- (The LL model appropriate for the low 
energy sector follows by effectively using a local inter- 
action, V{x,D) = VoS{x).) By virtue of a Bogoliubov 
transformation, H = Hq + Hy can easily be diagonal- 
ized. |-|With bosonic operators b'^\ the phase field ip{x) 
readsQ 



(p{x) 



q^Q 



.9(g) 
2L\q\ 



1/2 



exp{-iqx) sgn(g) [6j + b^g] , 



(8) 



with the g-dependent interaction parameter [where g = 
g{q = 2n/L)] 



g{q) = [l + e^Viq, D)/ttvf]-'^^ = g{-q) 



We then arrive at 



H ujqh\bg + pN 



(9) 



(10) 



with the plasmon dispersion relation LUg — vp\q\/ g{q). 
In terms of the chiral (right- or left-moving) phase fields 

{p = R/L = ±), 



(11) 



^pix) = [pipix)+^x)]/VA7T , 



the bosonized electron operato]p.|at a; = is ^'(0, t) cx 
Ep=±exp[27ri</)p(0,t)], implyingSH 

Goit) cx t-(9+9"')/2 e-iMt 

at long times. Hence one obtains the well-known expo- 
nent a = [g + g-'^ - 2)/2 governing the bulk TDOS. 

Let us now look at the static image potential experi- 
enced by an electron with charge — e held fixed at position 
a; = and < y < D due to its interaction with the QW 
electrons, Hj = —efdx n{x)V{x,y) = —e(j){y). Using 
Eq. (P), we get the fluctuating field 



cj,{y,t)^J2^i(y) (^«(*) 

q^O 

with couplings Xq{y) = -e[g{q)lq 
we shift the bosonic operators 
whence 



b-q{t)) 



(12) 



/27rL]i/2U(g,2/). Next 

Bq ^ bq- eXq{y)/ujq, 



2 



H 



qBlBg 



(13) which can easily be solved. For y = D, we finally get 



Here the static image potential 



V^miy) 



describes the energy gained by the plasmons relaxing to 
their equilibrium state in the presence of the additional 
electron. For the unscreened interaction (Q), one has 
g{q) = [1 + ^i^od^la)]"^^^ with the dimensionlcss pa- 
rameter ^ = j-KVp. Then the image potential (]lj) for 
{D — y) ^ a reads 



Vimiy) 



{D ~ y){H{D - y)/a] + 1/^} ' 



Therefore the static image potential is only logarithmi- 
cally suppressed against the bare Coulomb interaction, 
and hence is very strong. 

Next we turn to dynamical effects due to tunneling. 
We envision the latter as penetration through a rectan- 
gular barrier of width D. If the barrier is sufficiently 
thick and its transparency low, the main contribution to 
the tunnel current comes from electrons with momenta 
perpendicular to the QW. Therefore we effectively obtain 
a ID Schrodinger equation for the underbarrier motion 
■0(j/,t) of the tunneling electron (0 < y < D), 

tdt^Piy, t) = [-{2m)-'dl +ii' + U^ e^{y, t)]i^(y, t) , 

(15) 

where U is the work function of the tip/lead. In the 
absence of (jjly^t), the solution for an electron at energy 
^' reads ip{y,t) oc e~*^ where v — y/2U/m is 

an effective velocity related to the traversal time D/v. 
Under the WKB approximation, the dominant effect of 
the potential can be incorporated as additional phase 
factor ,Ej 



-ifi t — mvy — i9{y.t) 



(16) 



Linearizing the resulting WKB equation gives for 

\Hy,t)\ 



dt9 + ivdyO = —e(t){y, t) 



(17) 



supplemented by the boundary condition 0(0, t) = 0. 
This equation can be solved separately for each bosonic 
mode using the ansatz 

^(2/' ^) = E ^i)b-q{t) + w{y, OJq)bl{t)\ . 

<7#0 

From Eq. ( p7| ) and bq{t) = 6ge~*"«*, we obtain w{y,u!q) = 
w{y,-ujq) and 

{vdy - ujq)w{y,ujg) = ieXqiy) , 



0{t) = E H-^qHit) + W{u;q)b,q{t)] , (18) 



(14) with 



w{ujq) = (ie/v) dy Xq{y)ex.p[ujq{D - y)/v] . (19) 
Jo 

Because of the associated dynamic image potential, 
the electron wave function acquires the phase factor 
exp[—i9(t)] during the tunneling process. 

This effect can be properly incorporated by a modifi- 
cation of the tunneling Hamiltonian, 



Ht = r*^(0)e-*^co -l-h.c. 



(20) 



The Greens function determining the effective TDOS is 
now given by 

G{t) = (exp[i6't(i)]^'(0,i)^'t(0,0)exp[-i6'(0)]) . (21) 

Putting G{t) = Go{t)K{t), the contribution of the dy- 
namic image potential then gives rise to the factor 



where we introduce the functions 

C2(i) = 27r [{Mt)9{0)) + {9^{t)M0)) 



(22) 



where C2 is independent of p 
averages yields 



±. Doing the Gaussian 



Ci{t) = - '^[w{ujq)w{-ujq) + w^{-ujq) exp(-iwgt)] , (23) 



:[1 - cxp(-icj,i)]. (24) 



Since we are interested in the power-law exponent gov- 
erning the TDOS, we focus on the time-dependent parts 
of Eqs. ( p3| ) and (|2j), and do not explicitly compute the 
prefactor. For the interaction (|l|), numerical calculation 
of K{t) gives the result shown in Fig. |^ which is well ap- 
proximated by \K{t)\ = 1 + A{t) cos(rit) for long times, 
with oscillation frequency f2. The amplitude decays ac- 
cording to A{t) (X with (i ~ 1.15. This result is insen- 
sitive to the precise parameter values taken for ^,D/a, 
and v/vp. As a consequence, the power-law exponent 
a of the TDOS for small energy e remains unchanged 
by the dynamic image potential. Hence one can indeed 
obtain the LL parameter g from a measurement of a. 

This finding can be inferred analytically for a class of 
separable interaction potentials of the form 
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V{x,y) = Vo5{x)f{y) 



(25) 



where /(y) is an arbitrary function with ,f{D) = 1. In 
this case, we get V{q,y) = Vof{y), leading to g{q) = g. 
The time-dependent parts of Eqs. (p^) and (p4|) read 



D 



H 2 



dy f{y)e 



-vp{D~y)q/gv 



The asymptotic long-time behavior of Ci.2(i) can be ac- 
curately calculated in stationary-phase approximation. 
We find that Ci{t) decays faster than 1/t, while C2{t) oc 
1/t. Therefore, from Eq. (|2|), the TDOS exponent for 
small e remains unchanged. We expect this result to be 
correct and generic for arbitrary physically relevant in- 
teraction potentials. 

We conclude by summarizing our results. We have pre- 
sented a simple theory of electron tunneling from a tip 
or a lead into a strongly correlated ID metal, explicitly 
incorporating the finite traversal time and the dynamic 
response of the correlated metal to the incoming electron. 
We have solved this problem within a WKB-type approx- 
imation for different interaction potentials. Despite the 
presence of a strong static image potential, the power- 
law exponent entering the tunneling density of states is 
not affected by these effects, but completely determined 
by the correlation strength in the ID metal. 

We acknowledge financial support by the Deutsche 
Forschungsgemeinschaft under the Gerhard-Hess pro- 
gram. 
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FIG. 1. Tunneling from an STM tip into a Luttinger liquid 
(schematic). 
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FIG. 2. Time- dependent correction factor |-/^(t)| for the 
unscreened interaction (^. The parameter values are cho- 
sen as ^ = 1.72, D/a = 3, and v/vf ~ 1.48. These values 
should be appropriate for nanotubes, where 17 « 4 eV and 
vf ~ 8 X lO"" m/sec. The inset shows the decrease of the 
oscillation amplitude using double-logarithmic scales. The 
straight line is a guide to the eye only. 
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